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Abstract 

Perturbation of a propagating crack with a straight edge is solved using the 
method of matched asymptotic expansions (MAE). This provides a simplified anal- 
ysis in which the inner and outer solutions are governed by distinct mechanics. The 
inner solution contains the explicit perturbation and is governed by a quasi-static 
equation. The outer solution determines the radiation of energy away from the tip, 
and requires solving dynamic equations in the unperturbed configuration. The outer 
and inner expansions are matched via the small parameter e = L/l defined by the 
disparate length scales: the crack perturbation length L and the outer length scale 
I associated with the loading. The method is illustrated for a scalar crack model 
and then applied to the elastodynamic mode I problem. The crack front wave dis- 
persion relation is found by requiring that the energy release rate is unaltered under 
perturbation and dispersive properties of the crack front wave speed are described 
for the first time. The example problems considered demonstrate the potential of 
MAE for moving boundary value problems with multiple scales. 

1 Introduction 

Dynamic perturbation of a steadily travelling crack in a linear elastic medium is of fun- 
damental interest in fracture mechanics. The possible existence of edge supported mod- 
ulations in an otherwise straight edge raises questions about the stability of a steadily 
moving crack front. These waves, called crack front waves, have been the subject of intense 
scrutiny since they were first observed by Rice et al. pj and Perrin and Rice [Ij through 
numerical simulations of interactions of dynamic crack fronts with inhomogeneities. A 
theoretical framework for crack front waves was provided by Ramanathan and Fisher Q| 
using dynamic weight functions derived earlier by Willis and Movchan Ramanathan 
and Fisher [1] showed that a mode I (opening) disturbance can propagate along the crack 
front with a speed that is a function of the crack velocity and less than the Rayleigh wave 
velocity. Further numerical work confirmed the earlier findings [1[ and showed the explicit 
form of the (non-dispersive) crack front wave speed as a function of the crack speed v 
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The theory of Willis, Movchan and Ramanathan has been extended to other config- 
urations and experimentalists have sought evidence of crack front waves. Woolfries and 
Willis jsl gained insight by examining a scalar crack model. Woolfries et. al [|| gener- 
alised the dynamic weight function method to cracks in viscoelastic materials which was 
used by Willis and Movchan [?} to show that crack front waves decay in the presence of 
viscoelasticity. These theoretical studies considered perturbations of the crack edge in 
the plane of the moving crack. Obrezanova et. al 1(| and h| examined out-of-plane 



perturbations of a 2-dimensional crack, and Movchan et. al [12| considered the dynamic 
stability of a crack in a strip. Predictions of crack front waves have also been investigated 
through experiments. Crack front waves have been proposed as the cause of crack surface 
roughening in brittle materials [Hj]. Sharon et. al 13, EH, EH claim to have evidence 
of crack front waves in several experiments. Their conclusions are at odds with those of 
Bonamy and Ravi-Chandar [TtJ who used ultrasonic shear waves to distort dynamically 
growing cracks. They found that the perturbation of the crack front is a linear function of 
the wave amplitude, but the crack perturbation does not persist after the exciting wave 
has passed. Many questions on crack front waves remain to be answered. 

The purpose of this paper is to present a new, and, as we will argue, simpler, method 
for analysing the problem of a perturbed dynamic crack. The basic idea is that the size 
of the crack front disturbance is small compared with a macroscopic length scale, which 
we choose here to be defined by the dynamic loading. These yield inner and outer length 
scales and a small parameter e which is the ratio of the length scales. Our approach uses 
the method of matched asymptotic expansions (MAE) to split the elastodynamic problem 
into inner and outer sub-problems each of which is simpler than the entire problem and 
contains the physics appropriate to the region. As we will demonstrate, the inner problem 
is quasi-static and depends explicitly on the assumed form of the crack perturbation. The 
outer problem does not consider the fine scale of the perturbation directly, although it 
determines the radiation of energy from the inner region. The inner and outer solutions 
are related to one another through standard matching arguments. While the use of 
asymptotic methods in fracture mechanics is certainly not new, e.g. 18j . the method of 



MAE has not been used for dealing with complex moving boundary value problems of the 
type considered here. MAE is both useful and productive for this class of problem as it 
naturally splits the problem mechanically and mathematically. 

The layout of the paper is as follows. We begin in section [2] with a model problem 
demonstrating the matched asymptotic approach to studying the perturbation of a travel- 
ling crack edge. The steady solution and the scaling for the perturbation are introduced in 
section and the details of the MAE analysis are given in section [3j The same approach 
is then applied in section H] to the physically realistic case of a mode I crack travelling 
in an elastic material. The scaling, matching procedure and total solution for the elas- 
todynamic case are developed in a manner analogous to the scalar problem. The crack 
front wave dispersion relation is obtained and new results are presented for the dispersive 
behaviour of the crack front wave speed, with numerical examples given in section [5] 



2 The scalar model problem of a travelling crack 

In this section a model elastic problem of a dynamic perturbation to a travelling crack 
front is solved using MAE. The application is similar to that for the elastic crack but less 
complex. We therefore present it in detail in order to illustrate the general procedure. 
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2.1 Non-dimensional parameters and scaling 

Consider a function <j)(x,y, z,t), motivated by the anti-plane displacement field in mode 
III elasticity and satisfying a scalar wave equation 



V^-c"^ = 0, (1) 

where V 2 is the Laplacian operator, the subscript t denotes partial differentiation in time, 
and c is the speed of acoustic waves. A crack occupies the region 

x < vt + Lf(z,t), y = 0, — oo < z < oo, (2) 

that is, it grows at a constant speed, v, except for a perturbation term Lf, where the 
length L is the magnitude of the deviation and / is a dimensionless 0(1) quantity. For 
purely out-of-plane motions (displacement in the z direction) the crack faces carry zero 
stress, which we model here [see LLj as 

<f)y(x, 0, z, t) — 0, x<vt + Lf, (3) 

where <j) y = d<fr/dy, etc. The crack therefore opens at a rate that is approximately 
constant, and it is the deviation from the steady state that is of interest. The opening is 
assumed to be caused by a steadily translating symmetric loading on the crack faces, 

. P f x — vt\ 
<f> y (x,±0,z,t) = -p — — . (4) 



I r \ I 

Note that we have assumed that the non-dimensional forcing function p has argument 
that scales on a length scale I, which is taken to be much larger than the perturbation 
scale L. (For simplicity we interpret this to also imply that p(x) is zero in the vicinity of 
the crack tip.) We define the small parameter e as the ratio of the two length scales: 

e = L/l. (5) 

It is convenient to move to a coordinate system fixed in the (constant) moving reference 
frame, and to non-dimensionalise the initial boundary value problem. Thus, we define 
dimensionless independent variables as 

x-vt y z ct 

— -*> y-^ T -*> y-*> (6) 

together with 

and the problem reduces to analysing the system 

a 2 4>xx + 4>yy + 4>zz - 4>tt + 2f (f) xt = 0, (8a) 

(f) y (x,±0,z,t) = p(x), x<ef(z,t), (8b) 

where 

a = s/l - v 2 . (9) 

The loading function p(x) is, by definition, zero close to the crack tip, and we assume that 
the crack growth speed is subsonic so that a remains real. The system (jSj) still needs to 
be supplemented by a crack growth criterion, discussed below. We emphasise that here 
and henceforth all parameters and variables are non-dimensional. 
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2.2 Constant running crack 

It will prove useful later to first consider the solution to (JSJ) for zero perturbation in 
the crack tip position, f(z,t) = 0. We start with a constant running point loading 
applied symmetrically to the crack faces at distance unity (i.e. at distance I in the original 
coordinate system) behind the crack tip: 

<j> y (x,±0,z)=6(x + l), x<0. (10) 

Here 5(x) is the Dirac delta function. Note that, with a symmetric loading, will be zero 
ahead of the crack tip, and so this is a standard two-part boundary value problem which 
may be solved in a variety of ways; a solution of (jHJ) and ffTO]) is therefore 



(j) y (x,y,z) = Re< — — — — - -. \ . (11) 

[ ir(x + lay + Ijyx + lay J 

It may be checked, using the identity (x — iO)^ 1 — (x + zO)" 1 = 2iri5(x), that ( TTTT) indeed 
satisfies (TTU]) . Expanding about the crack tip, and integrating, the solution for is found 

as 

{—1 °° (— l) n i] 

Note that in (fTl]) and fTT^j) . Re and Im denote the real and imaginary parts of the ex- 
pressions to their right, and the curly brackets are henceforth omitted to maintain clarity. 
More generally, we consider a symmetric loading along the crack faces such that 



'y 



[x,±0,z)=p(x), x<0 (13) 



with, again, the assumption that p(x) is zero in a region close to the crack tip. Based on 
the line load solution, the near-tip field is 



{x, y, z, t) = Im P {x + icty) * [l + — (x + lay) + ...], (14) 



where 



— CXD — OO 

The amplitude Pq defines the strength of the square root 'displacement' behaviour behind 
the moving tip and Kq is the analogous scalar 'stress intensity factor', such that the stress 
ahead of the tip is 

(f> y {x,0,z,t) = (1 +mx + ...) . (16) 

V27TX 

Note that m = — 1 for the line load limit (flOl) . and the subscript in P an d K Q denotes 
the displacement and stress coefficients for the unperturbed crack tip at x = 0. Finally, 
the scalar energy release rate is 

G = P K (17) 

where, in general, 

P = lim \(-xy 1/2 (P(x,0,z,t)} , K =]3mU/x<l>y(x,0,z,t)]. (18) 
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3 Asymptotic solution; scalar model 

In the previous section the length scale of the forcing was taken to be much longer than 
that of the disturbance to the crack tip. This disparity introduces the small parameter 
e into the problem (in f!8bl) ) which can be usefully employed to solve the problem using 
the method of matched asymptotic expansions The non-dimensionalisation carried 
out previously has cast the initial boundary value problem into the 'outer-field' form so 
we shall have to derive an 'inner' coordinate system too. In the two coordinate systems, 
asymptotic expansions will be deduced, containing sequences of unknown coefficients, and 
these will be determined by matching together the respective terms. 

3.1 The outer expansion 

We begin with the following ansatz for the outer expansion 

(f>(x, y, z, t) = (o) + e0 (1) + e 2 (2) + . . . , (19) 

where the superscripts in brackets here and henceforth refer to the expansion function at 
the order in e indicated. Each term in the expansion must satisfy eq. (JHal) as well as the 
crack face conditions 

0( n) (x,±O,^,t) =5 n0 p(x), x<0, n>0, (20) 

where 5 n o = 1 if n = and 5 n o = otherwise. The leading order term is the solution for 
the unperturbed crack, 

0(°) = ImP o [s 1/2 + ^ 3/2 + -..], (21) 
where the complex variable s is 

s = x + iay. (22) 

The complex form of the solution automatically satisfies a homogeneous crack boundary 
condition close to the tip. Note that, in terms of the outer coordinates x and y the crack 
perturbation is very small, hence the reason for locating the centre of the coordinate 
system at the origin. We will return to the next term in the outer expansion, (f)^ 1 ', once 
we have considered the leading terms in the inner expansion. 

3.2 The inner expansion 

We introduce dimensionless inner variables X and Y, 

x = eX, y = eY, (23) 

and let Q(X, Y, z, t) = 0(x, y, z, t). The wave equation expressed in the inner variables is 

« 2 $xi + $yy + 2ev<S>xt + e 2 ($ zz - $ u ) = 0, (24) 

and the edge is now located at X = f(z, t), Y = 0. The inner ansatz is motivated by the 
scaling and by the form of the leading order outer solution: 

$ = e l/2 $(1/2) + e 3/2 $(3/2) + _ _ _ _ (25) 
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The leading order term is easily shown to take the form 

$(V2) =i mA ^S^ 2 , (26) 

where the inner complex variable S is 

S = X - f + iaY , (27) 

is centred on the perturbed crack tip in order to capture the correct singular behaviour at 
the shifted edge. The potential <l>( 1//2 ) satisfies the governing equation ([2"4"]1 with e set to 
zero (i.e. a scaled form of Laplace's equation) and its derivative is zero on the crack faces 
as required. The unknown A^ 1 ' 2 ^ is found by matching the inner and outer expansions. 
More specifically, we rewrite the inner expansion in terms of the outer variables. Thus, 
taking $ up to (^(e 1 / 2 ), expanding it in terms of the outer variables up to O(e ), and 
using S = e _1 s — /, we obtain 

= lmA^s 1 / 2 . (28) 

Similarly, the expansion of <fi up to O(e ) when expressed in terms of the inner variables 
up to C(ea) is 

0{o,|} = £ i/2 Im p ^ x + ia y)l/2 _ (29) 

Comparing (|28|) and (|29|) . we see that they are equivalent if A^- 1 ' 2 ^ = P$. 

The next term in the inner expansion satisfies, according to (1241) and (|25|) . 

a^S + *?f = -2i>4>r • (30) 
The solution is a sum of the particular integral and a general solution. The former is read- 
ily found using 3>jJt = \Poft im S~ 3 ' 2 , combined with the identity (a 2 d 2 x + dy)Sg(S) = 
4a 2 g'(S), where S is the conjugate of S. Adding the appropriate general solutions, we 
determine 

$(3/2) = Im p o r A (3/2 )>S 3/2 + 5(8/2)5,1/2 + V g S -l/21 ^ 

L 4cr 
where A^ 3 ' 2 \ £?( 3 / 2 ) are as yet unknown. The expansion of $ up to 0(e 3 ' 2 ) in outer 
variables up to O(e ) is therefore, 

${|> > = \mP [s l l 2 + A { - 3 ' 2 h 3 ' 2 ]. (32) 
Similarly, from (|T9|) and (l2Tj) . 



1 >-ImP [ S 1 / 2 + ^ 3 / 2 ]. (33) 



Comparison of the latter two expansions, which should be identical by the matching rule, 
implies A^ 3 ' 2 ^ = m/3. 

To summarise, the inner solution has been determined to 0(e 3 ' 2 ) as follows, 

$ = ImPo k 1/2 S^ 2 + e 3 l 2 {-S 3 ' 2 + B^S 1 ' 2 + —J t SS- 1 ' 2 )} + G(e 5 / 2 ) . (34) 
L v 3 4cr 

The single real- valued coefficient _g( 3 / 2 ) remains unknown, and will be found in the next 
section. Finally, we note for future reference, 



$ {|,i} = ImP 



(s 1/2 + js 3 ' 2 ) + e [ - \f s^ 2 + (B^ - ^fy' 2 + g rss^ 2 } 



(35) 
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3.3 Wiener-Hopf analysis 

In order to proceed to the next order in the outer solution, which is required if we are to 
match with the inner expansion from (134)) it is necessary to consider the perturbation in 
the wavenumber frequency domain. The most general form of f(z,t) can be constructed 
from the solution 

f(z, t) = Im f e i{ - kz - ut) = Im /„ e^"*"^ , (36) 

where k = k/u (assumed to be less than unity) is the edge wave slowness and fo is a 
constant. We seek possible solutions of the outer system of equations with no forcing but 
which display the singularity represented by the s -1 / 2 term in (1351) . We therefore assume 

<t>W{x,y,z,t) = 1mq(x,y)e i "(« e - t \ (37) 

where q(x,y) satisfies, according to (JBJ), 

ce q X x + q yy + (^ — k )q — 2iujvq x = 0, — oo < £ < oo, ?/ > 0, (38) 

and 

q y (x,0) = 0, x<0, (39a) 

q{x,0) = 0, x>0, (39b) 

lim (x 2 + y 2 ) 1/4 g(x, y) < oo. (39c) 

x 2 +y 2 — >0 

One type of general solution is found by a standard analysis involving the Wiener-Hopf 
method. Thus, let 

oo 

g(X ' y) = i / ^(0^^, (40) 

— oo 

where 7 follows from (l3~8l . 

7 (o = (e 2 + ^ 2 -(i-o 2 ) 1/2 . (41) 

The integration contour runs along a strip in the complex £ plane, T> say, which contains 
the real line except that it is indented above the point — A + and below the point A_, where 
these correspond to the branch points of 7. The Riemann surface of 7 is selected so that 
.Re 7 > 0. We further define 

7 = 7+7- , (42) 

with 

7±(0 = a (£ ± A ± ) 1/2 , A± = (1 - k V) 1/2 ±- 2 . (43) 

By definition, the function 7 + (7") is analytic in the upper (lower) half plane containing 
the common strip of regularity T>. Generally, the superscript notation ± indicates hence- 
forth functions analytic in these overlapping half-plane^- The boundary conditions on 
x < and x > imply, respectively, 

7 g(0=T+(0, (44a) 
q(0 = W~(0- (44b) 



1 However, the reader is reminded that the subscript on the constants A± refers only to the choice of 
sign in ([15])! 
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Thus 7 W (£) and T + /7 + are equal to one another and hence by analytic continuation 
and Liouville's theorem must be also equal to a constant, say q . Therefore, g(£) = 

?o/7"(0- 

The near tip behaviour of q(x, y) follows from the behaviour of the transform at large 
£, thus, as r —> 0, 



DC 



q(x,y) 



d£ e~ u ( i Z x+ ^\ a v) 



27ra J {t-iO) 1 / 



2 



l+~^ay(X+-X-) sgn£+m/(A + +A„) 2 ^ + 



(45) 

This can be evaluated using the identities 



oo 



1 f ^e-^K^+lcps) 2 n' e l7T ' i 

^ 7 (e _ i0) i/ 2e , [l,^e]=-^^(-O w lH^[Mn-|)g,«oos(n-j)g], 

— oo 

(46) 



for n > 0, where s = x + zcn/ = |s|e . Hence, 



g e l7r/4 / . ,_ 1/2 . . 1/2 . 6 

q[x,y) = ^= s ' sin — h «a>A_ s 7 sin - 



+ [^(K - X-)\s\ 1/2 - ^(A + + A_) 2 |s| 3 / 2 ] sin (9 cos t + . . . ) . (47) 
Thus, from (1371) . the second order term in the outer expansion is 

0(D = i m (\ s \-\ sin i+L\s\^ [(A++3A_) sin -+(A + -A_) sin -1 ^ e^"^ +0(|s|§ 

a^/^fuJ\ 2 4 LV ' 2 2 J / 

(48) 

We are now ready to complete the matching with the inner field of ((35D. This 

requires that the coefficients of |s| -1 / 2 sin |, {s] 1 ^ 2 sin | and Is] 1 / 2 sin y, are identically 
equal. The former implies that 

g = \a^e~™ / 4 P /o, (49) 
while the \s\ l l 2 sin | terms match if 

£(3/2) = imp + * ( A+ + 3 A _)l /oe i(fe- 1) 
2 8 

= Im lj + ^ (- 2 " * V ) ""I A*"*™ + (50) 

where ( l43]) 9 has been used to express it in a form that will be useful later. Finally, we note 
that the Is] 1 ' 2 sin y terms automatically match on account of the identity (A + — A_) = 
2v/a 2 . We are now ready to consider the complete MAE solution. 
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3.4 Solution and discussion 

The strengths of the perturbed singularities at the perturbed crack edge are defined by 

P = e" 1 / 2 lim iXl'^iX, 0, z, t)=P + eP (P (3/2) + -^f t ) , (51a) 
xio 4a 

K ee r^liml^^H^O.M) = K + eK (B^ - ^Lf t ), (51b) 

where the values follow from the leading order terms in the inner solution in (|34p using 
( 1501) . Consider the variations from the values for the steadily propagating crack: AP = 
P — Po, AK = K — K Q and AG = G — Go, where G = PK is the perturbed energy release 
rate, then the relative changes are 

AP vf t AK vf t AG 

= eq * t + e — -, = eq * t — e — -, = eq * t, 52 

P y J 2a 2 ' K y J 2a 2 ' G K J 

where * denotes convolution, and the transform of g(t) is 

g(u,k) = m + ia~ 2 (LU 2 -k 2 a 2 ) 1/2 . (53) 

The expressions for AK and AG are identical to analogous ones derived by Woolfries and 
Willis (1999) (eqs. (2.5), (2.7) and (2.18) of their paper). They also agree with prior work 
by Rice et al. (1994) for the special case of m — 0. 

Crack front waves are possible if the phase speed c p = uo/k and the crack speed v 
together lie inside the sonic cone, i.e. they satisfy 



c? t =v 2 + c? p <l. (54) 



In that case the transform g becomes 



c 



1-v 2 

Zeros are possible only if m is positive and satisfies 



g(u,k)=m-\k\— 5-. (55) 



< m < \k\/y/l ~ v 2 , (56) 



in which case g(u, k) has a unique zero at uo = a\Jk 2 — a 2 m 2 . This solution has been 
discussed before, but we note briefly some properties in the wavenumber domain. The 
phase speed as a function of k is monotonically increasing from zero at the cutoff or lowest 
wavenumber possible, k = myl — v 2 , to its asymptote c p = y/l — v 2 as k — > 00. The 
crack edge waves are therefore dispersive, and one can, formally at least, define a group 
speed (velocity), c g = duj/dk. This is related to the phase speed by 

CpCg = 1-v 2 , (57) 
and hence, for finite frequency /wavenumber 



v 2 + c 2 g > 1. 



(58) 
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This type of dispersion is anomalous in that the group speed is normally associated with 
the energy propagation speed. However, there is no strict connection, and in fact, the 
notion of energy propagation for edge waves has not been defined. As a result, we may 
define and calculate c g although its physical significance is not entirely clear. 

It is also of interest to examine the final solution for the near field, as it provides 
the complete behaviour in the neighbourhood of the moving crack front. Equation (IMj) 
implies that the inner solution can be expressed 

$ = P„([l + e (\g . / + ^)]p" 2 sm | + _ si „ 36) + 0(f3/2) _ (M) 

where p, G are polar coordinates relative to the perturbed moving tip: x — ef+iay = pe . 
We note the appearance of the p 1//2 sin ^ term, which affects the stress singularity ahead 
of the crack, but has a different angular dependence from the standard near tip field 
p 1 / 2 sin ® for a steadily moving crack. 



3.5 Properties of the MAE scheme and simplification 

Before considering the elastic crack in detail, we note some general features of the matched 
asymptotic analysis for the scalar problem. The sequence of terms derived was of 
(fit]) -> §W of $25)) -> $( 3 / 2 ) of (El]) -> of (@BD. The latter was derived as an 
eigensolution using the Wiener-Hopf method, and was necessary in order to complete the 
matching of the single remaining coefficient in $( 3 / 2 ). Note that the particular solution 
of the inner term for the forcing in eq. ( 130]) was not required for matching. In fact, the 
terms in AP and AK that involve ft cancel in the final expression for AG, i.e. g of ( 153]) . 

Thus, we calculated the perturbed values of P and K separately, although the quantity 
of interest, G, which is their product, turns out to be simpler. This suggests a more direct 
procedure, which follows by noting that G has alternative expressions: 



G = PK = ^1 (1 - V 2 )- 1 ' 2 K 2 = y^- (1 — V 2 ) 1 ' 2 P 2 , 



(60) 



where V = v + ef t is the velocity of the crack edge. The identities (160]) are a consequence 
of the general form of (fl5j) i for speed V, 

K=^-{l-V 2 f/ 2 P (61) 

As a consequence the perturbed energy release rate can be expressed in different ways, 

AG AP AK AK v „ AP v „ 

G P K K a 2 P a 2 

Only one or other of P or K needs to be considered in order to calculate the energy release 
rate. In practice, it is simpler to compute P as it is the field variable calculated in the 
MAE procedure. Restricting attention to the crack face X < 0, Y = +0, and using eqs. 
flMD and (|35J>, we have 

$(-|A|, +0, z, t) = P [eV 2 (l + ep^Xl 1 ' 2 - e 3 / 2 ^|X| 3 / 2 ] + G(e 5 / 2 ), (63) 
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and 



\X\,+0,z,t) = P 



\x 



1/2 



-\ x m 



X 



^ 1/2 + (pi-f/)N 1/2 ] 



(64) 



where 



4cr in 



(65) 



At the same time, the outer expansion evaluated on the crack face is, from (1331) . fl37|) and 

ED, 



D iu)(i\x\ 



A. 



m = ^M ft -j\^ + '^^^ J ^^{i + ^ + CK^)]. (66) 



The integral can be evaluated using (I46p . Matching the 0(e) terms in (1641) and (I6"6"|) . the 
| a? | 1//2 singularity gives the identity (1491) for q , while the next term yields 



Pi 



/ 



(67) 



The perturbed energy release rate then follows from eqs. (j62l) . ( |65l) and ( j67l) . 

In summary, we have presented two methods to complete the matching - the first and 
more general applies the matching to the field <fi(x, y, z, t) for all x and y near the origin, 
and the second only matches the field on the crack faces, X < 0, Y = 0. The two methods 
are equivalent because the final stage of the matching requires only a single real-valued 
coefficient, and therefore matching of the field along a line amounts to matching over 
an area (in x and y). However, £?( 3 / 2 ) is not the most relevant quantity for the purpose 
of calculating the perturbation in energy release rate, although this coefficient can be 
found using the second approach, from eq. fl65l) once the coefficient pi is known. These 
same general features are repeated in the elastic case. In particular the calculation of the 
perturbed energy release rate will be achieved directly once the matching is completed. 



4 The elastic crack 

We now turn to the more realistic problem of in-plane perturbations of a steadily propa- 
gating mode I crack in an isotropic elastic solid. It is important to distinguish between 
two distinct types of 'crack edge waves'. The first is the analogue of the Rayleigh surface 
wave, that is, an infinitesimal wave of particle displacement confined to the edge and 
decaying away from the edge. The second is a wave-like modification of the edge itself. 



Achenbach and Gautesen [20] demonstrated the nonexistence of the first type of edge 
wave on a stationary crack edge. It can be shown using their arguments combined with 
the analysis in the Appendix, that this result extends to to a steady propagating crack. 
That is, there there are no localised solutions of this form for any crack velocity that is 
subsonic relative to the Rayleigh wave velocity. 
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4.1 Scaling and asymptotic expansions 

The unperturbed crack lies in the plane y = 0, with infinite edge moving steadily at 
speed v, located at x = 0, — oo < z < oo, in the convected coordinate system. We 
adopt the same scaled coordinates as in the scalar problem, with the distinction that 
the normalisation in time is with respect to the Rayleigh wave speed cr, so that the 
non-dimensional variables are 

x-vt y z c R \ v 

-,-,—t}^(x,y,z,t), >v. (68) 



I I I I J Cr 

The inner and outer expansions for the displacement field u = (u x , u y , u z ) T are 

u(x, y, z, t) = u (0) + eu (1) + e 2 u (2) + . . . , (69a) 
U(X, Y, z, t) = eVSud/^ + e 3/ 2u (3/2) + _ ? (69b) 

where the inner variables X and V are as the same as before, eq. ff23l) . The expansions in 
fl69|) are the analogs of those in eqs. ffl9l) and ff25l) for the scalar problem. The boundary 
conditions on the crack faces are that the traction vector cr = (a xy , a yy , u zy ) T vanishes, 
except for the loading, defined below. 

We introduce the non-dimensional speeds associated with the two bulk wave speeds 
in an isotropic elastic medium 

vj = cj/cn, I = L,T, (70) 

where cl = y(A~+ 2/i)/p and ct = \J ^/ po are the longitudinal and transverse wave 
speeds, A and /i are the Lame moduli and po is the mass density. The displacement field 
is represented by three potentials: 

u = V0 L + V A <p T e 3 + V A V A (pTH^s, (71) 

where the L and T are associated with the longitudinal and transverse waves, respectively, 
and TH is transverse horizontal. The second transverse potential, <Pthi is zero in the 
absence of the perturbation. The effect of changing to the moving coordinate x is that 
each 0/ satisfies a modified wave equation 

1 v 

afyl^x + 4>I,yy + 4>I,zz - ^2 + 2 ^2 = °> I = L,T,TH (72) 

with 



vj vj 



a I = a I (v) = yjl-v 2 /vl I = L,T, (73) 

and axH = ®t, vth = vt- The asymptotic analysis will be performed in terms of the 
potentials, represented as a 3-vector for both the outer and inner regions: 

(74) 

\$thJ 

respectively, and which are expanded as 

4>{x, y, 2, t) = (o) + e0 (1) + eV (2) + • • • , (75a) 

F, Z, t) = e V2$(l/2) + ,3/2^(3/2) + _ _ _ _ (?5b) 
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4.2 Steadily propagating crack 

We first consider solutions for the steadily propagating crack in an isotropic elastic 
medium. The unperturbed solution has 4>th = and the potentials 0l and 0^ may 
be represented by analytic functions of two complex variables 



sj = x + ia>iy, 



L, T. 



(76) 



We introduce the real valued potentials, i/> n , which we express as functions of the complex 
variables sl, St, 



ip n (s L ,s T ) = 



I2V2ttD(u + §)(n+ \) 



Re 



i2ai s T 




(77) 



where 



D(v) = 4a L ax 



a 



2 \2 
TJ i 



(78) 

and D(l) = based on the normalisation of the wave speeds with respect to the Rayleigh 
speed. The associated vectors of displacement and traction are 



1 



Re 



-1 



ItxD 



Re 



a T )s L 



-2icti 




i2otL Sy ' * 



[1 + 4) s n L 2 ( 1 + af, I + i2a L s T 




(79a) 



(79b) 



We note in particular, that the zero traction conditions on the crack face are satisfied. 

These fundamental elements may be used to describe a steadily propagating mode I 
crack, with zero applied shear on the crack faces, i.e. a xy (x, ±0) = cr zy (x, ±0) = for x < 
0. For instance, the solution for a pair of travelling line loads with <J yy (x, ±0) = —5(x + 1) 
is 



. 1 
a yy (x,y) = — Re 



7T 



D 



[l + s T )s^ /2 



'1 + a 



2 \2 



sl s 



1/2 



(80) 



The near tip fields, in this case, can be expressed terms of the basis functions explicitly 

as 

5>irY n , tr= J-X)(-1) B S„. (81) 



u 



7T 14 — ' V 7T 

n=0 n=0 



More generally, the loading 

a yy (x, ±0) = —p(x), x < 0, (82) 
implies a near tip stress expansion identical to that for the scalar problem (see (TLB"]) ) 



(f> — K Q (?/> + mi/?! + . . . ) cr TO (x, 0) = (i _|_ mx 



v/2 



7TX 



+ ...) 



S3) 
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Here Kq is the stress intensity factor, H is the Heaviside (step) function, and 

o o 

K °=)J- J dsi-Sy^pis), m =-±-yjl J ds(-s)^ 2 P (s) 

— oo — oo 

This solution serves as the basis for the perturbed crack. 



54) 



4.3 Perturbed crack: inner and outer analysis 

The leading order term in the outer expansion (j 75 aft is, from fl83l) . 



= K [V> 0l, s t ) + mip^SL, s T ) + 0(4 , s 



3/2 3/2 

"t 



55) 



Turning to the the inner expansion we note that the potentials satisfy modified wave 
equations 



1 



a 



]®i,xx + $/,yy + e2-$ 7 , X t + e 2 (<S> IjZZ - =0, I = L, T, TH 



The leading term in the inner expansion (175bj) follows from fl85l) as 

<S>W = K *I> (S L ,S T ), 

where the inner variables are 

5/ = X - f + , I = L,T. 



(86) 



(87) 



The potentials defined by ( 187|) satisfy the governing equations ( 186|) . and the traction 
conditions on the crack faces. 

The potentials in the next term of the inner expansion satisfy the inhomogeneous wave 
equations 

(89) 



^2^(3/2) .(3/2) _ _«i!_ffi(l/2) r r T T rr 



These may be solved in the same manner as before, as a sum of homogeneous solutions 
plus the particular solution, 

^ = K [aW^S^St) + B^ (S L ,S T ) +vM,(S L ,Sr)], (90) 
where the final term is the particular solution given by 

/ (i+<4) 5 q1/2\ 



*J>*(S L ,S T ) 



1 



fiV27iD 



Re 



V 



a L v L 
Hag 



2iotT c G^/2 







(91) 



/ 



Matching the outer term (p^ '^ with $^2'°i implies that A^/ 2 ^ = m. 
The complete inner solution is therefore 

$ = K [eV 2 ^ (S L ,£V)+e 3/2 (mt^ 

(92) 
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from which we can express the expansion of the inner in terms of the outer variable as 
^{f' 1 } = K (^ (s L ,s T ) +rm/>i(si,s T ) + e\^f%l)_ x (s L ,8 T ) 

+ (5(3/2) _ i m f)4 ( aL ,8T)+vft1> m (8 L ,8T)]\ ■ (93) 

These two equations are the equivalents of eqs. (|34|) and fl35l) for the perturbed scalar 
crack model. We note that, as in the scalar problem, the issue is reduced to finding a 
single real constant, B^ 3 ^ 2 \ 

4.4 Matching and final solution 

Based on the experience with the scalar problem, we assume that the 0(e) outer solution 
is time-harmonic and of the form 

0« = Re e M^-t) u_ I d £ e -i^x ^(i) ^ ^ 

where 

k = 1/cp (95) 

denotes the slowness for an assumed edge disturbance with phase speed c p . The crack 

edge perturbation itself is defined by (|36|) . The transform </> (£, y) is derived in the 
Appendix, to within a multiplicative factor. The latter is determined by the near-tip 
expansion, which is expressed via the large |£| behaviour of the transform. Thus, from 
the Appendix, we have 



<fp) =R ee luj{KZ - t) — / d£ 
2vr 



£3/2 



f-i(l + a\ )e-^l Q ^ (sgnO (l + 0{Z~ X ))\ 
-2a L e-^ aTy (l + C(r 1 )) 

o(r 2 ) 



= Re e M«-*) 6 ^=- jj,D [1>_ X (8 L , s T ) + 0(sf, sf )] , (96) 

where the identities ( 1461) have been used. The coefficient 60 is found by matching the 
coefficient of with the 0(e) term in ( 1931) . to give 

e-"' 4 A'„/o 

"» = _ yl~~^D~' <97) 



Hence, the near tip expansion is 



The next order term can be evaluated using the explicit expressions in the Appendix, and 
in principle the near tip expansion can be continued. However, the leading order term 
suffices for our present needs. 



Perturbation of a moving crack edge 



1(3 



The energy release rate for the perturbed crack front is 

G = PK = FK 2 = F~ l P 2 . 

Here, P and K are defined 



(99) 



P = e~ 1/2 lim \\X\- 1/2 U y (X,0,z,t 



K = e 1/2 lim [\X\ l/2 E yy (X, 0, z, t)] , (100) 

X [0 



where U y and H yy are the leading order term in the inner expansion of displacement, 
(I69bh . and stress, and |2l| 

with V — v + ef t . Based on the lessons learned from the scalar problem, we use the final 
expression in ( 199]) to find the perturbed energy release rate as 



AG AP F'(v) 
G ~ P F(v) U 



(102) 



Following the procedure for the scalar problem, we can obtain AP by matching the 
crack opening displacement. Thus, 



U y (-\X\, +0, z, t) = -?L [e l ' 2 {\ + e Px )\X\ 1 / 2 - e^\Xf 2 ] + 0(e 



5/2 > 



'2ll 



and 



u}*' 1} {-\X\,+0,z,t) = 



2tt 



l/2_!^| x |3/2 



)+e [i/lxl-^+^-^/JIxl 1 / 2 ] 



where now 



The outer expansion for the crack opening is 



(103) 



(104) 



(105) 



u. 



{1,1} ^0 



.r 



^-||x|^ +£ Ree«--)^ / ^^(i + l+^r"))], (106) 



where the complex number a is defined by the expansion of u y 1 ^ in eqs. ( 1A.16I) and ( 1A.17I) . 
Matching the the |a;| _1 / 2 singularity in the O(e) terms in (I104p and (I106[) gives the identity 
(I97P . while the \x\ 1 ^ 2 coefficient yields 



m 



Pl - -f = Im iuaf e i <- kg - ut \ 



(107) 



Hence, the transform for the perturbation in energy release rate, defined by (1521) ^. is 

F'(v) 



g(tu, k) = m + iiola + itu 



F(v) 



(101 
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This is the fundamental dispersion relation, which will be discussed further in the next 
section. We note that a similar expression has been obtained by 0, eq. (4.10)], and by 
(author?) 0, eq. (3.21)0. However, their expressions differ from (11081) in the sign of the 
second and third terms. This distinction is important for computing the dispersive form 
of the crack front wave (see section (151)), but is immaterial in the non-dispersive limit of 
m/\k\ — > 0. Finally, we note that the real parameter _E?( 3 / 2 ) follows from eqs. ( 11051) and 
(11071) . and the complete inner expansion to 0(e 3 ^ 2 ) is then given by eq. fl93|) . 




A 




0.2 0.4 v 0.6 0.8 1 



Figure 1: The limits of the parameter A of eq. (1112j) as a function of the crack speed. The 
Poisson's ratio of the material ranges from 0.45 to —0.95 in increments of 0.1. 



5 Numerical results and discussion 

We discuss solutions to g(u, k) = for < v < 1, where v is the normalised crack speed 
relative to the Rayleigh wave speed. Crack front waves with phase speed c p = uo/k lie 
inside the Rayleigh sonic cone if 

c 2 = c 2 + v 2 < 1, (109) 

where c is the total speed of the disturbance in the fixed stationary frame. It may be 
shown [ll] that the transform g of (11081) is real valued for edge wave speeds inside the sonic 
cone. Furthermore, using the change of variables of Ramanathan and Fisher (3|, also see 
[(|, we have 

\k\-^(u k) — A - 2^HZ + /dp [2^-cV+ P )M P ) mo) 



2 The parameter m used here is the same as that of (author?) [9|] but twice the value of the quantity 
ui m (author?) 0). 
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where 

Pip) = tan ( 2 ), (HI) 

and 

A=|fc|- 1 m=— . (112) 

27T 

Here A is the wavelength along the crack, and A is therefore the non-dimensional ratio of 
the wavelength to the length scale of the dynamic loading. At a given crack speed, crack 
front waves exist for a range of this parameter: A min < A < A max , where the lower and 
upper limits correspond to c = 1 and c — v, respectively. Thus, 



AmlnM = - ylzlpi _ / m^p-^-pWp) (113a) 



T 



. dp v 2 (3{p) . . 

Amax(f) = 7j = " - - / —7- ov, /2 ^ - (H3b) 



Figure [TJ show A min and A max as a function of the crack speed for the range of elastic 
materials, characterised by the Poisson's ratio v = (\v 2 L — v^)/{v\ — Vj<), which has the 
permissible range —1 < v < 0.5. We note that A min is always negative and decreases 
monotonically to a finite value as v — > 1. Note also that A max increases monotonically 
from unity at v = and goes as A max « 2/ \Jl — v 2 as v — > 1. This behaviour is to be 
compared with the scalar problem (see eq. (1551) ) for which A mill = and A max = l/y/l — v 2 . 
The dependence upon Poisson's ratio is not strong, and therefore for the remainder we 
consider the single case v = 1/3. 

The phase speed is plotted in Figure [2] as a function of the crack speed for a range 
of the non-dimensional wavelength A. This extends the results of Morrissey and Rice [y] 
who presented the A = curve for the total speed c, and of Willis and Movchan [3] who 
computed the A = curve for the phase speed. Figure [3] shows the total velocity c as a 
function of v for the same set of A. The total speed is always close to unity for negative 
values of A. Thus, 0.995 < c < 1, 0.987 < c < 1 and 0.961 < c < 1 for A = -0.4, -0.2 
and 0, respectively. Figure 0] plots the phase speed as a function of A for two values of 
the crack speed. The range of A is consistent with Figure [TJ and note that the lower and 
upper limiting values of the phase speed correspond to c = v and c = 1, respectively. 

By analogy with the scalar problem, the group speed of the crack front wave travelling 
along the edge is defined c g = du/dk. It may be shown to be 

,2„,2 , ^2„ 



A/ 2 v L f dp (cV + c 2 p - 2p 2 )(3(p)\ 



9 P c p \(l-v 2 )VT^ (v 2 -v 2 )^l^ J ir {p-v 2 )\p-c 2 Yl 2 ^p) ' 

(114) 

and is plotted in Figure [5] for a range of values of the wavelength-loading parameter A. 
Note that values of c g in excess of unity occur, which is highly anomalous. However, in 
contrast with the scalar case, the group speed does not always satisfy the inequality (j58|) . 
although it holds more often than not . 
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Figure 2: The crack front phase speed c p for v = 1/3 vs. the crack propagation speed for 
seven values of the non-dimensional wavelength: A = —0.4, —0.2, 0, 0.2, 0.4, 0.6 and 0.8. 




Figure 3: The total speed of the crack front wave c = ^Jv 2 + c 2 vs. the crack propagation 
speed for seven values of wavelength: A = —0.4, —0.2, 0, 0.2, 0.4, 0.6 and 0.8. 
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Figure 4: The dependence of the crack front phase speed c p on the non-dimensional 
wavelength A at two values of the crack speed, v = 0.6 and v = 0.8. The upper limits of 
c p in each case are a/1 — v 2 . 



2 




Figure 5: The group speed c g of the crack front wave vs. v for the same set of values for 
A as in Figure [31 
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6 Conclusions 

We have shown how to solve the scalar model of a perturbed travelling crack edge by using 
the method of matched asymptotic expansions. The general procedure developed for the 
model elastic problem was then applied to the elastodynamic mode I crack. The MAE 
solution employs outer and inner expansions which are matched via the small parameter 
e = L/l defined by the disparate length scales: the size L of the edge perturbation and the 
length scale I associated with the loading. The key to the solution is a homogeneous outer 
eigensolution which displays unphysical singularities at the unperturbed crack tip. The 
reason for this is that the crack is shifted, although the unphysical eigensolution is not 
directly applicable at the crack edge. No unphysical singularities are present in the inner 
solution, or in the complete matched solution. We have used the MAE scheme to derive 
the dispersion relation for the crack front wave speed, and have shown its dependence on a 
range of parameters, including crack speed, wavelength, and elastic properties. Once the 
MAE solution is found to second order in both the inner and outer fields, the dispersion 
relation is obtained by requiring that the energy release rate, which depends on the inner 
solution, is unaltered. 

The intent of this paper has been to show by example that MAE provides a natural 
methodology for solving problems of this nature. Future work will develop the procedure 
to situations that are not easily solved by other means, for instance the perturbation of a 
dynamic crack front with a cohesive zone at the travelling crack tip. In this case the size 
of the cohesive zone defines the inner length scale. 

Appendix 

A Solution of the O(e) outer problem 

We consider the symmetric boundary value problem defined by the boundary conditions 
on y = 0: 

(A.la) 
(A.lb) 

®xy = ®zy = 0, — OO < X < OO. (A.lc) 

This defines a mixed boundary value problems on the split regions x < and x > 0. The 
solution is assumed in the form ()94p . and henceforth we will omit the term e luJ ( KZ ~ t ) an d 

concentrate on the transformed quantities u^(£, y) and <x^ (£,?/). 

The symmetry of the problem implies that we need only consider the half space y > 0. 
The transform of the potentials in this half-space may be represented by 



a yy 


= o, 


x < 0, 


Uy 


= 0, 


x > 0, 


— &zy 


= 0, 


— OO < X < OO 



(iu^Ce-^y 



= (zcj)- 1 I £ e -™ ( A .2) 



and hence the displacement is 
u« (£, y) = A(-£, iry L , n) T e~™ + [B(i lT , £, 0) T + C(-£k, i 1t ^ 1% - £ 2 f] e~™, (A.3) 
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where 



il = (e +k 2 - \(i - vif) i/ \ 7T = (e +k 2 - - o 2 ) 1/2 . (a.4) 

v L v T 

The square root functions are taken so that Re7^ > 0, Re7T > 0, thus ensuring decay in 
the half-space y > 0. 

The associated stresses follow from eqs. (172!) and (1A.4|) . The transformed traction, 
when evaluated on the plane y = 0, becomes 



(iw/O-^C&o) 



-2i lL i -e-i 2 T 



-27j,ac 



/A 

B 



(A.5) 



It is also useful to list the analogous relation for the on-plane displacement transform, 
from (1A.3|) with y — 0, 



u«(£,0) 



«7i £ «Tt« 
k 7t " £ 2 



D 



(A.6) 



We are now ready to consider the boundary conditions. 

The condition that the stresses a xy and o zy vanish on the plane y = implies that 

&xy 0) — and a% ■ (£, 0) = 0. The first and third rows of the vector relation (1A.5I) 
then allows us to eliminate any two of (A, B,C) in favour of the remaining quantity. 

Substituting the result into the equations for a yy (£, 0) and tt^(£, 0), we deduce that 



4 1} (£,o), 



7l(1-<) 2 " 
where -R is the (modified) Rayleigh function 



(A.7) 



(A.8) 



The boundary conditions (lA.la|) and (lA.lbl) imply that a yy \^, 0) is a (+) function, i.e. 
analytic in the upper half of the complex £-plane, and My (£, 0) is a (— ) function. Define 



Q(0 = Q + (0<T (0 



(A.9) 



where D is defined in ( 1751) and £ = — ?7 + and £ = r/_ are the zeros of the modified Rayleigh 



function i.e. the roots of £ 



;i-o 2 = o, 



77 ± = -^-(v/l-K 2 (l-i; 2 )±t;). 
1 — t? z ' 



(A.10) 
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Thus, 

Q(0-1, Kl-oo, (A.ll) 

and the functions Q ± (0 are defined unambiguously by requiring that they both have 
this property. Using the above mentioned analytic properties of the transforms we may 
rewrite (IA.7I) as 

where, by the usual analytic continuation arguments, is an entire function of the com- 
plex variable £ and 

7^(0 = a L (£ ± A i± ) 1/2 , A L± = - «^|c* 2 ) 1/2 ± . (A.13) 

Assume that the stress and displacement behave as 

0) = ^((-x)- 1 / 2 ), 0) = 0((-*r 3 / 2 ), (A.14) 

near the edge. This implies that the transforms behave as £ -1 / 2 and f 1 / 2 as f tends to 
infinity, and therefore the entire function E is in fact a constant. We assume the constant 
is non-zero otherwise <fi^ = and the matching cannot be accomplished. 

We note that the only non-trivial solution to ( 1A.7I) has an unphysical singularity near 



the edge of the crack. This removes the possibility of infinitesimal localised disturbances, 
analogous to Rayleigh waves on a free surface. Achenbach and Gautesen [20| considered 
the possible existence of travelling waves on the edge of a stationary crack, and showed that 
neither symmetric nor anti-symmetric solutions can exist. Equation flA.7j) generalises their 
result to steady propagating cracks, and shows the non-existence of symmetric localized 
modes. A similar analysis can be done to show that anti-symmetric modes cannot exist 
on the propagating crack. 

In general, the near tip behaviour follows from the expansion of the potentials as 
|f | -> oo. Using 

oo 

■If.. In Q(Q \ 



we have 



Q-«) = exp(-y dC^j, (A.15) 



4 1) «,0) = ^[1 + | + O(r 2 )], (A.16) 



where, from flA.12j) and flA.13j) 



oo 

I A --i/ dchl « (c) - (A - 17) 



This is used in section (jlj) to obtain the perturbed energy release rate. The detailed form 
of the inner expansion and the near tip field can be obtained by noting that the explicit 
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form of the potentials are 



B 

w 

which imply 



c 



4 1} (e,o) 



f(e+i 2 T )(e-K 2 -i 2 T )+2e^\ 
-^7L(e 2 + ^ 2 -7^) 

\ 2k7 L 7 T / 



B = -2za L C 1 {l + ^ + ---)u ( y\t,0), 



_ai) -] + ...)4 1) (e,o), 



a-. 



r "-lfHi + i[i+3] + ...)«? , (f.o)- 



as 



Equations (IA.16I) and (lA.19j) can then be used to derive the near-tip fields. 



(A.18) 



(A. 19a) 
(A. 19b) 
(A.19c) 
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